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which is a necessary condition to apply the time-sharing argument to 
an AVC subject to a state constraint. 

The fact that (1) does not imply (16) also appears to prevent one 
from applying Ahlswede's elimination technique [1], [7], to AVC's 
subject to a state constraint; thus, our result in Theorem 1 does not 
appear to be sufficient to allow us to prove the conjecture that the 
region in (2) is exactly the capacity region when state constraints are 
present, either by the elimination technique or by time-sharing. 

Corrections to [6] 

The preceding observations require us to make the following 
corrections to our prior paper [6], 

In [6, Theorem 5.8] the words 41 closed convex hull,** should be 
replaced by "closure." 
m In [6, Section V-C] we can no longer conclude that equation (5.8) 
is correct. Instead, all that we can conclude, using notation defined 
^ in [6], is that the closure of 

®l/ 2 (p** Q*, W a ) U ^ /2 (/7*, <?*, W a ) 

is a subset of C{W a , 1/2) and that C(\V a , 1/2) is a subset of the 
right-hand side of (5.8). 

Acknowledgment 

The author is grateful to an anonymous reviewer of an earlier 
version of this correspondence for reminding him that the time-shar- 
ing argument can not be applied because (1) does not imply (16). 

Appendix 

Lemma; Let p be any distribution with p(x) > 0 for all xeG. 
For any qe ^(G), H(p*q) - H(q) = 0, if and only if q is the 
uniform distribution. 

Proof: Define a distribution on G X G by setting P XY (x, 
y)=P(x)q(y - x). Then 

l(XA Y) = H( Y) - H(Y\ X) 
= H(p*q) - H{q). 

Clearly, H(p*q) - H(q) = 0, if and only if X and Y are inde- 
pendent. Since p(x) > 0 for all x, X and Y are independent, if 
and only if q(y - x) = q(y) for all x, y e G. Thus, X and Y are 
independent, if and only if q is uniform. □ 

Theorem A: Let W be a general multiple-access AVC. If 
p e ^(X) and q e ^(Y) are such that (recall paragraph 2 of Section 
III) 

I(XaZ) > l(SAZ), for all re^(S), (A.l) 

and 

I(YAZ\ X) > I(SaZ\ A r ), for all re9 L (S), (A.2) 

then every pair (/?,, R 2 ) satisfying 

0<R X <I L {XAZ) and 0 < R 2 < I L { Y A Z | X) (A.3) 

belongs to the deterministic-code average-probability-of-error ca- 
pacity region under state constraint L. 

Proof; This theorem can be proved by making trivial modifi- 
cations to the proof of [5, Theorem 5.5]. A similar observation was 
made in [6, Section V-A], though it was not pointed out there that in 
this case the modifications do not require that the channel be 
nonsyrnmetrizable [6, Definitions 3.3, 3.5, and 3.7]. D 
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Abstract — Constellation shaping is extended to provide shaping gains 
without resorting to high-dimensional constellations. This is accom- 
plished by dividing the constellation into unequal sized constellations, 
and selecting these constellations on an equiprobable basis. A design 
example is provided, demonstrating the simplicity and power of the 
approach. 

Index Terms— Coding, modulation, shaping, signal constellations, 
and nonequiprobable signaling. 



I. Introduction 

Coding schemes for transmission of data over the Gaussian 
channel have been used over the last decade that lead to improved 
performance [l]-[4]. The most popular approach is to use coset 
codes, and to attain high code rates coupled with good performance, 
the shift has been to use higher dimensional constellations. One 
result of moving to higher dimensions is the ability to achieve what 
has been called shaping gains. This is due to the reduction in 
average symbol energy that can accompany the use of a constella- 
tion whose boundary is not an A/-cube. In particular, as the constel- 
lation becomes more spherical, it enforces a nonequiprobable distri- 
bution on signal points drawn from a constituent two-dimensional 
constellation. It has been shown [5] that in the limit as N — ► oo, an 
A^-sphere can achieve 1 .53 dB of shaping gain, and will enforce a 
truncated Gaussian distribution on the constituent 2-D constellation. 
Attention has been focused on these shaping gains, as they can be 
achieved independently of any gain due to the use of a coset code. 

Multidimensional constellations with significant shaping gain were 
first described by Conway and Sloane in [6], Other constellations 
with shaping gain and simple decoding methods were described by 
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Calderbank and Sloane in [7]. Forney [8] extended the method in 
[7] to more general lattice partitions. Subsequently, Forney has 
applied this approach to constellations constructed from trellis codes 
in [9], and Eyuboglu and Forney have applied it to nonideal 
channels [10]. 

Recently, it has been proposed that the most direct route to 
achieve shaping gains is to use lower dimensional constellations and 
enforce a nonequiprobable distribution on the signal points [11]. 
Using this method, Calderbank and Ozarow have been able to equal 
the best shaping gains found using higher dimensional constella- 
tions, and yet have better peak-to-average power ratios (PAR) than 
found in [8]. 

The method of Calderbank and Ozarow divides the constellation 
into T subconstellations of equal size. A shaping code is then 
applied to select subconstellations according to a given probability 
distribution. Signal points within subconstellations are equiprobable. 

In this correspondence, we consider extensions of the work of 
Calderbank and Ozarow. In Section II, we describe shaping as 
proposed by Calderbank and Ozarow, and derive an alternate method 
of shaping that offers higher shaping gains. We also extend the PAR 
calculations to this new approach. This method is essentially the 
same as that first described in [2], but found through a different 
route. In Section HI, we generalize the results of Section II, and 
describe our shaping method as the dual of source coding using 
prefix codes. In Section IV, we compare shaping using prefix codes 
to previously published shaping gains, and demonstrate that higher 
shaping gains are possible with potentially little complexity. In 
Section V, we describe the major drawback of this method, namely 
that a random number of bits are now transmitted per channel use. 
We describe two methods proposed to alleviate the problems associ- 
ated with this fact ([12], [13]), and offer an application where it may 
have little impact. Section VI provides an example of signaling 
using nonequal subconstellation sizes, and demonstrates that for 
very little complexity shaping gains of 1 dB can be achieved. 
Section VII summarizes the results. 

II. Shaping Gains 

Shaping gains refer to the reduction in average signal power by 
constellation design and signaling choice. In [11], the total shaping 
gain is contributed to by two separable factors. The first is the 
shaping gain of the constellation chosen. This gain varies between 0 
dB (for AAcube) to 1.53 dB for a spherical constellation in an 
infinite number of dimensions. The second factor contributing to the 
shaping gain is what Calderbank and Ozarow call the biasing gain, 
and comes about through the use of a nonequiprobable signaling 
scheme. 

A, Equal Sized Subconstellations 

In [11], the nonequiprobable signaling is achieved by first divid- 
ing the constellation into T subconstellations, all of equal size. This 
subdivision is defined by first defining a fundamental region 0t . The 
construction continues by scaling the fundamental region T times, 
to form a constellation as in Fig. 1. Thus, the T copies of the basic 
region are a, 01 , * - • , a r _ l 0t, With the nested sequence of 
scaled regions defined in this manner, we may define the total 
constellation, Q, as the union of T subconstellations, Q f , i — 
0,---,r- 1, where fl,= Q Pi (a, + 1 ^/a, -^). Signal points in 
each subconstellation are used equiprobably, and subconstellation fl, 
is used with frequency f - r 

In this correspondence, Calderbank and Ozarow make use of the 
continuous approximation to compute average power. This approxi- 
mation assumes that the constellation consists of a set of continuous 




Fig. 1- Decomposition of a rectangular constellation into subconstellations. 

points. This gives results that are easily calculated, and yet are 
accurate for large constellations. Using this method, the biasing gain 
normalized to two dimensions is given by [11], 



TatC/o , * * * » fr- i ) — 



2(2/A/)/^( J r 0 ."-./r-i) 



£#--o[('+0 



2/N+ I 



;2/N+ 
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where T is the number of subconstellations, N is the number of 
dimensions, and f t is the probability that subconstellation /, / = 
0, — ,T — 1, will be chosen. The reader is referred to [11] for 
details. 

In [11], the increase in the PAR due to constellation expansion 
and the use of a shaping code was found to be 



PAR, - 



(2) 



where T and the f { are as in (1). Thus the peak-to-average power 
ratio of the nonequiprobable signaling constellation is given by 



PAR Q « PAR 2 x PAR, 



(3) 



where PAR a is the peak-to-average power ratio of the fundamental 
region, 01, and PAR n is the peak-to-average power ratio of the 
total constellation, Q. 

Calderbank and Ozarow are able to achieve their shaping gains 
without increasing the PAR beyond that found using higher dimen- 
sional constellations. In fact, they find that their method yields equal 
shaping gains, but with better PAR figures. In some cases, the 
improvement is by a factor of 2 or more over the other methods. 

B. Nonequal Subconstellations 

We are left to wonder if this construction is optimal, namely, the 
splitting of the region into T equal sized subconstellations. [11] 
argues that equal sized subconstellations are better than subconstel- 
lations with diminishing size, i.e., with K(fl 0 ) > K(fi,) > • • - > 
^(ft r -i)» where K(Q,) represents the volume of region 12 ( -. They 
provide an example with V(£l;) = 2K(Q /+1 ), where this is indeed 
the case. However, other relations between the subconstellations are 
possible. 

We begin by generalizing (1) to find the biasing gain for subcon- 
stellations that are not equal. We first note that in [11], the average 



BNSDOCID: <XP 2871 44A > 



8097 IEEE Transactions on Information Theory 
38(1992) July, No. 4, New York, US 



IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 38, NO. 4. JULY 1992 



1349 



transmitted signal power (normalized to two dimensions) is derived 
from the expression: 



1 



(/+ 0 fa 



£/ y p 0 = ('+0 1/A >o. 



(4) 



where P Q is the average power of the fundamental constellation Q 0 , 
Recall also that the number of points in each subconstellation is 
equal . 

Now let the subconstellations be constructed as follows. Again we 
scale the fundamental region T times, yielding T copies of the basic 
region ® , a, @ , * ■ ■ , ct T _ , ^ as before. However, now let the 
number of points in each subconstellation, Q h / = 0,'-'J- 1, 
vary. In fact, let us define 



(5) 



The quantity (kj - k J _ l )/k T _ i then represents the fraction of 
signal points in the yth subconstellation, and the following holds 
true: 

K T~\ 7=0 

Solving for tj yields 



kj - kj, i 



Thus, the average power in the yth annular subconstellation is 



tjPo = 



kj kj_ , 



(?) 



(8) 



Given we choose subconstellation Q { with frequency f h then the 
total average power is 



P = 



T- 1 /-2/A/+I _ 

E\> K j-\ f 

j = 0 K j K J~l 



(9) 



Selecting subconstellation j increases the rate by log 2 (£ y - — 
kj_ ,), thus, using nonequiprobable signaling with unequal subcon- 
stellations yields a total rate of 



E // log (*/ - *,_ ,) + H(f 0 , • • • , />_,) 



(10) 



bits per channel use. Equiprobable signaling at this rate yields an 
average power of 



2 VN 



T- I 



£ /, log (*, - *,._,) + H(f 0 r • * , f T - ,) 



00 



This yields a biasing gain normalized to two dimensions of 
7A/(/o'*'*>/r-i) 

2^ N [T.T^f t log (*/-*/-,) +ir(/o,-"./r-i)] 



k 2/N+\ _ 



(12) 



Continuing with the approach of using different sized subconstel- 



lations. we find that the increase in the PAR is given by 

k 2 f N 
k r- i 



PAR, = 



r r-i 



k{ — k t _ , 



(13) 



where T, the f h and kj are as in (12). 

At this point we pause to note that using the new construction, it 
is theoretically possible to achieve slightly higher shaping gains and 
slightly better PAR'S than using equal subconstellations. However, 
as will be seen in Section IV, the advantage of the new construction 
is to be found when using realizable shaping codes. 

III. Generalization of Construction 

We have been motivated so far by extending the method of 
Calderbank and Ozarow to nonequal subconstellation sizes. Through 
this approach, we have arrived at an expression that yields the 
shaping gain for a wide range of relationships between constella- 
tions. A more general description can be found, and in fact points 
out the relationship between this method and that described in [2]. 

In [2], the idea of shaping gains using prefix codes was described. 
This method can be explained as follows. If we desire to send an 
average of ft bits per two-dimensions, then start .with a two-dimen- 
sional constellation with M > 2 & points, and assign a /3,-bit label to 
the yth point such that 



(14) 



where P- is the probability of selecting the yth point. If our 
selection of labels is done in such a way as to satisfy the Kraft 
inequality, 



E2-»/Sl, 

J 



(15) 



then there exists a prefix code of alphabet size 2 with the integers 
/3,, 0 2 > ' * " 2 s code-word lengths. If the signal points are selected 
by applying the prefix code to the incoming bit stream, and if the 
bits are equally likely, then the probability of selecting signal point 
j is now 2~^. If we assign longer labels to points farther from the 
origin, then a Gaussian distribution can be approximated. A simple 
way to implement this scheme would be to group points together 
(i.e., define a subconstellation as described in the previous section) 
and assign fixed-length labels to points in the same subconstellation. 
This description of shaping may be regarded as the dual to the use 
of prefix codes in source coding. 

This formulation can be seen to be more general than that 
described in Section II. The simple implementation previously de- 
scribed may also be seen to include as a natural progression the case 
where kj - & y _ i decreases with increasing y, as with the general- 
ized multidimensional cross constellations described by Forney and 
Wei ([5]), to kj - kj_ x = K y a constant, as in Calderbanlc and 
Ozarow, to where kj - kj_ , increases with increasing y\ as de- 
scribed in Section II. In fact, one may generalize the multidimen- 
sional cross constellation construction even further by splitting the 
constellation into more than two subconstellations. The construction 
then proceeds along the same lines as found in [5], and preserves 
the important property of transmitting an equal number of bits per 
multidimensional signal point. However, (12) may be seen as a 
simple way to calculate the biasing gain provided by the general 
approach if we assign fixed-length labels to points in the same 
subconstellation. 
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IV. Implementation of Coding 

A. Shaping Codes 

Calderbank and Ozarow propose the use of a nonlinear code to 
provide the shaping. The code is defined as follows. Let B n ^ m be 
the set of all binary n-tuples v with Hamming weight wt(u) < m. 
Then, the rate of the code is given by 



Shaping Bits 



(16) 



This code can be used to select two subconstellations with proba- 
bility f 0 and by selecting constellation Q 0 when a zero is 
present in the codeword, and selecting constellation Q l when a one 
is present. The probability / 0 is given by 



(17) 



Note also that when using a code, the term H(f Qi /,,' • * , fr- 1) m 
(1) and (12) is replaced by R(B n ^ m ), the rate of the code. 

Using these codes for T = 2, Calderbank and Ozarow are able to 
achieve practical biasing gains of 0.51 dB, using a B 205 nonlinear 
code. However, encoding such a code can be difficult. For more 
practical length codes, where a lookup table can be used, gains are 
less than half a dB. Table I lists some of these codes for T = 2 and 
T = 4, with associated code parameters and shape gains. Note that 
for T = 4, [1 1] combines two B n m codes to select between the 
four subconstellations. 

When we turn to nonequal subconstellations, we find that the best 
B n m code for any T is a B n „ code. That is, we choose subconstel- 
lations equiprobably. This corresponds to the use of a rate R = 
log 2 (7*) code. If we were using equal sized subconstellations, then 
a rate R — iog 2 (T) code simply implements an equiprobable sig- 
naling scheme. However, with nonequal subconstellations, using 
equiprobable subconstellation selection still yields nonequiprobable 
signaling on the whole. 

The shaping code now reduces to labeling the T subconstellations 
with binary n-tuples of length log 2 (7*). We use Iog 2 (7") bits at the 
input of the coder to select the subconstellation, and signals are then 
selected from the subconstellations on an equiprobable basis. Given 
that the number of symbols in the yth subconstellation is M Jt then 
the total number of bits required to select any symbol in subconstel- 
lation j is log 2 (Mj) + log 2 (T) bits. This fits in well with the 
description of shaping as the dual to source coding using a prefix 
code, as given in Section m. Here, P { - log 2 (Mj) + \o$ 2 (T), if 
signal point i is in the yth subconstellation. 

Since there are a different number of symbols in each subconstel- 
lation, a serial to parallel converter can be used to provide the bits 
required to select the symbol from the subconstellation. In Fig. 2, 
we have illustrated a block diagram that implements a shaping code. 

Using T = 2, and splitting the two constellations into nonequal 
sizes, we find that the optimum ratio of the outer constellation size 
to the inner constellation size for equiprobable subconstellation 
selection is 3 to 1 (i.e., k l = 4, k Q = 1). This yields a biasing gain 
of 2/ VT (0.62 dB), a PAR 2 of 1.32, and a constellation expansion 
ratio of 1.15. (As an aside, the biasing gain is exactly the ratio of 
the density of the lattice A 2 to the density of lattice Z 2 ([7]).) Thus 
we can theoretically achieve a shaping gain of 0.82 dB at a cost of a 
15 per-cent increase in constellation size, and a 32 per-cent increase 
in PAR. All of this is achievable without significant increase in 
complexity. Table II lists relevant parameters for T = 2 and T = 4 
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Fig. 2. Block diagram of shaping code implementation. 



TABLE I 



N 


T 


n 


m x 


m 2 


72 


CER 2 


PAR 2 


7 


PAR 


2 


2 


6 


2 




0.35 


1.19 


1.29 


0.55 


2.58 ; 


2 


2 


12 


3 




0.45 


1.24 


1.38 


0.65 


2.76 


2 


2 


18 


4 




0.51 


1.26 


1.41 


0.71 


2.82 


2 


4 


6 


1 


3 


0.43 


1.55 


1.71 


0.63 


3.42 


2 


4 


12 


2 


5 


0.60 


1.50 


1.73 


0.80 


3.46 * 


2 


4 


18 


3 


7 


0.68 


1.48 


1.73 


0.88 


3.46 



equiprobable selection of unequal subconstellations. Note that these 
in all cases have better biasing gain and PAR 2 than the codes listed 
in Table I for the same T. In Table m, we list the shaping gains 
available using multidimensional equiprobable signaling for various 
lattices as described in [11]. The Leech lattice is denoted A 24 , y 5 
denotes the shaping gain of the Voronoi region of the lattice, and 
PAR 2 denotes the peak-to-average power ratio of the constituent 
2-D constellation. Note that the equiprobable subconstellation ap- 
proach also outperforms each of these lattices. 

B. Realizable Shaping Code Constellations 

One drawback to the use of non-equal subconstellations is the 
requirement (in binary signaling) that the number of signal points in 
each subconstellation must equal a power of two. Under this restric- 
tion, the k t in (2) can take on only a discrete set of values. (This 
restriction can be relaxed somewhat at the expense of implementa- 
tion complexity). 

In Fig. 3, we have plotted the realizable biasing gain for T = 
2, 4, 8, 16 versus the increase in PAR over the constituent constella- 
tion, PAR 2 . For comparison, in Fig. 4 we have plotted biasing gain 
vs. PAR 2 for T= 2,4 for both approaches. Note that the new 
method in all cases has both higher biasing gain and lower PAR 2 . In 
both figures, (U) refers to the unequal subconstellation construction, 
and (E) refers to the equal subconstellation construction. In Fig. 5, 
we have plotted the maximum shape gain for both methods and 
several lattices vs. the peak-to-average power ratio. The parameters 
for the lattices are found in [11] and [8]. We assume a spherical 2-D 
constellation for the shaping codes. The shaping gain includes both 
the biasing gain and the constellation shaping gain. The PAR is the 
PAR of the constellation for methods (U) and (E), and is the PAR of 
the constituent 2-D constellation for the lattices. The reader is 
referred to [8] for an explanation of the notation for the lattices. 
Here, we see that the new method achieves higher gains at lower 
PAR than any of the other approaches. We also point out that the 
PAR of the lattices can be reduced by using peak constraints. 

V. Encoding and Decoding Issues 

Here, we briefly describe decoding of the shape codes. In [1 1], a 
level slicing approach is proposed, where the received signal is used 
to decide which subconstellation was in use. If a B n m code is used, 
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TABLE II 

Biasing and Shaping Gains for Equiproejable Selection of Subconstellations for T ~ 2 and T = 4 
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r 
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A*3 


72 
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1.15 
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TABLE III 

Shaping Gain Parameters of Some Popular Lattices 



Biasing Gain as o Function of PAR2 
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Biasing Gain 


as a Function of PAR2 




o 



(U) T= 2 
(U) T= 4 
(U) T= S 
(U) T=16 



1.5 



2.0 



2.5 



PAR2 



Fig. 3. Realizable biasing gains for new construction and T = 2, 4, 8, 16. 

then after n symbols have been received and decoded by subconstel- 
lation, we can form the codeword used to provide shaping. This 
codeword isthen used to find the actual data, by using table lookup, 
or some other block decoder method. 

Since this method relies on the use of a nonlinear code, and the 
minimum distance of the code is one, this approach can lead to 
increased errors. For example, using T — 2, and a B 20 ,5 code, if 
one symbol is received in error causing us to select the wrong 
subconstellation, then with high probability we will decode incor- 
rectly. This incorrect decoding can affect up to 14 bits at the output 
of the decoder. However, an important property of this approach is 
that each symbol transmitted represents the same number of bits, 
thus limiting their effect to a finite length block. 




Fig. 4. Comparison of realizable biasing gains for both constructions. 



However, using the new method, since the subconstellations each 
represent a different number of bits, errors on reception can lead to 
bit insertions and bit deletions. This is a significant problem, and 
must be addressed in any practical system. Here, we describe two 
proposed methods to deal with this, and propose an application 
where insertions and deletions do not cause unsolvable system 
problems. 

The first approach we describe was proposed by Gallager in [12]. 
This document describes a system that implements shaping through 
the use of a secondary channel. The constellation is divided into 
subconstellations, similar to that as described in Section II. How- 
ever, the extra bits needed to select points in the outer subconstella- 
tions are drawn from a secondary asynchronous source. The trans- 
mit sequence is designed in such a way that if bit deletions or 
insertions occur, the receiver automatically compensates by insert- 
ing or deleting, respectively, the same number of bits into the 
received asynchronous data stream. This ability to resynchronize, 
while not avoiding the erroneous reception of data, makes the 
approach attractive from a realizable system viewpoint. 

The ^synchronization is accomplished by use of a transmit 
buffer, and a specific schedule for the minimum number of bits to be 
held in the buffer at the end of each symbol transmission. If, after 



BNSDOCID: <XP 2871 44A > 



8097 IEEE Transactions on Information Theory 
38(1992) July, No. 4, New York, US 



IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 38, NO- 4, JULY 1992 



the end of a modulation period, there are not enough bits in the 
transmit buffer to meet the scheduled minimum, then the shortfall is 
made up from the secondary asynchronous data source. At the 
receiver, if an error is made such that it believes it has received, 
say, one more bit than was actually transmitted, then, because of the 
scheduled minimum number of bits in the transmit buffer, it will 
believe the extra bit came from the secondary data source, and send 
it to the secondary data sink, thus re-acquiring synchronization with 
the primary data source. A similar process occurs if the receiver 
believes one less bit has been transmitted than actually sent. In this 
case, the transmitter will at some time transmit an extra bit from the 
secondary data source, but the receiver will send it to the primary 
data sink, thus again resynchronizing the primary data. For a full 
description of this ingenious approach, the reader is referred to [12]. 

This approach points out one of the problems that variable-rate 
signaling faces, that of buffer management. Since a random number 
of bits are transmitted during each interval, some buffering must 
eventually be employed. During any transmission interval, the 
encoder may request either more or fewer bits than requested in the 
previous transmission interval. How this data is provided without 
the input buffer overflowing or becoming empty is of serious 
concern. Any system using variable-rate signaling will require some 
scheme to provide input data to the encoder in a timely fashion. 

The second method proposed to combat this problem, is to 
partition the input data into constant length blocks, and then force 
the system to transmit this information using the same number of 
symbols each time. This method is described by Calderbank and 
Klimesh in [13], and uses what are known as balanced codes. 
Balanced codes have the property that every codeword has an equal 
number of zeros and ones. Using a simple construction provided by 
Knuth in [14], Calderbank and Klimesh are able to construct rate 
2 p /(2 p + p) balanced codes with a particularly simple method of 
encoding and decoding. Using this approach on nonequal subcon- 
stellations, with codes of high rate, they are able to preserve almost 
all of their shaping gain. For example, it is possible to construct a 
rate 256/264 balanced code. Since the number of zeros and ones are 
equal in each word, we have f Q = /, = 1/2. Using this code as the 
shaping code, with nonequal subconstellations having k Q = 1 and 
k x = 3, yields a biasing gain of 0.44 dB, only 0.09 dB less than the 
biasing gain offered using a rate 1 code. In this case, we would have 
to partition the input data into blocks corresponding to 264 symbols, 
with half of the symbols drawn from the inner region. Thus every 
block would transmit 1321og 2 ( | ft 0 | ) 132 x 2log 2 ( | Q Q |) + 
256 bits. 

A place where the shaping described in Section II could be 
applied most simply is in an ARQ system. Here, the input data is 
collected into a group of a fixed number of bits to form a packet. A 
cyclic redundancy check code is applied during transmission. Upon 
reception, the CRC is used to check for errors. If errors have 
occurred, then a retransmission is requested. Headers and trailers 
using fixed constellations would also simplify this system, allowing 
the receiver to detect the starting and ending sequences that delin- 
eate a packet. Again, the insertion or deletion of bits in this case is 
not a fatal error, since a retransmission can solve the problem. 
However, if the system is using synchronous multiplexing between 
ports, then bit insertions and deletions may prove to be an insur- 
mountable problem. 

VI. A Design Example 

We conclude this correspondence with a simple example of 
shaping using nonequal subconstellations. In this example, we select 
T = 4, with k Y = 3, k 2 = 7, and £ 3 = 15. We select as a funda- 
mental region $t the 16-QAM constellation, with the second, third, 
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Fig. 5. Comparison between shaping gains offered by several 
approaches. 
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and fourth subconstellations containing 32, 64, and 128 points, 
respectively. Thus, our total constellation will consist of the 240 
points illustrated in Fig. 6. Note that the subconstellations are 
outlined in this figure. 

We will use an equiprobable subconstellation selection shaping 
code. With four subconstellations, the rate is 2 bits /channel use. 
Thus the biasing gain is theoretically 0.87 dB. Calculating the 
average power of each subconstellation, we find that (assuming the 
points lie on an odd half-integer lattice) P Qo = 2.5, P Qi = 10.25, 
P u = 25.375, and P n ^ = 56.125. Thus, the total average power is 
P s = 23.6. Note that the constellation, as it stands, has a shaping 
gain of nearly 0.2 dB. 

To compare for gain purposes, we note that given the constella- 
tion in Fig. 6, the average number of bits transmitted per channel 
use is 7.5 bits. Computing the baseline average power as in [2], we 
find P u = 2 7 5 /6 = 30.17. Thus, the total shaping gain (including 
biasing gain and constellation shaping gain) is I0log 10 (30.17/23.6) 
= 1.07 dB, equaling the sum of the biasing gain and the shaping 
gain of the constellation. This is a significant gain for little complex- 
ity. If we were to use balanced shaping codes as proposed by 
Calderbank, then, using two rate 1024/1034 codes, the loss in 
biasing gain would be 2 20/l034 = 0.06 dB, yielding an overall gain 
of 1.01 dB. 

Similarly, PAR 2 = 1.76 and-CER = 1.33. A better comparison 
than PAR 2 is to use the ratio of the PAR,, of the shaped constella- 
tion to the PAR U of the unshaped constellation. This yields 

PAR. , , 

= 1.22. (18) 

PAR„ 
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Fig. 6. 240 point constellation used in the design example to achieve 1 dB 
of shaping gain. 

Thus, the increase in PAR is actually less than that indicated by 
PAR 2 . 

VII. Conclusion 

We have described an approach to achieving shape gains that, 
while not completely new, has largely been ignored. We have 
demonstrated the power of the approach in that it achieves higher 
shaping gains than several other schemes, and with less constellation 
expansion and peak-to-average power increase. However, there are 
serious system issues to be resolved before the method can be 
successfully applied. 
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Weight Distribution of Cosets of 2-Error-Correcting 
Binary BCH Codes of Length 15, 63, and 255 

Paul Camion, Bernard Courteau, and Andre Montpetit 

Abstract— The weight distributions of cosets are not known for the 
class of binary 2-error-correcting BCH codes of Length n = 2 m - 1, m 
even (the nonuniformly packed case). By using the graph theoretical 
concepts of combinatorial matrix of a code and /--partition design 
introduced in previous works, we obtain these weight distributions of 
cosets in the three particular cases of length 15, 63, and 255. We observe 
that in these three cases the number of distinct weight distributions is a 
constant equal to 8. 

Index Terms — BCH codes, weight distribution of cosets, automor- 
phism groups, coherent partitions, partition designs, regularity number. 

I. Introduction 

In [6], the notion of r-partition designs was introduced to extend 
to more general situations some properties of partial difference sets 
[3], [7] which is the particular case corresponding to 2-partition 
designs. If a linear code C admits a r-partition design, which means 
that the set of columns of a parity matrix of G is a class of the 
partition, then the theory developed in [6] gives the weights of the 
orthogonal code and allows us to effectively calculate the distance 
matrix of the code whose rows give the coset weight distributions of 
C. The aim of this paper is to use the methods of [6] to calculate the 
distance matrix in the three unknown cases of 2-error-correcting 
binary BCH codes of length 15, 63, and 255. 

The class of binary 2-error-correcting BCH codes of length 
n — 2 m — 1 and codimension 1m has been extensively studied. In 
1960, Gorenstein, Peterson, and Zierler [11] obtained that these 
codes are quasiperfect, the covering radius being 3. In 1967, 
Kasami [12] has determined the weight distributions of these codes. 
When m is odd, the orthogonal code has three nonzero weights 
whereas when m is even it has five nonzero weights. So when m is 
odd the binary 2-error-correcting BCH codes of length n = 2 m - 1 
are uniformly packed [15] by a theorem of Goethals and Van 
Tilborg [10], [16] (see also [1] and [2]) and the theory of Van 
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